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This paper starts with a systematic rederivation of the statistical thermodynamic equations of 
motion for dislocation-mediated plasticity proposed in 2010 by Langer, Bouchbinder and Lookman. 

It then uses that theory to explain the anomalous rate-hardening behavior reported in 1988 by 
Follansbee and Kocks, and to explore the relation between hardening rate and grain size reported 
in 1995 by Meyers et al. A central theme is the need for physics-based, nonequilibrium analyses in 
developing predictive theories of the strength of polycrystalline materials. 


I. INTRODUCTION 

In an earlier publication [I| (LBL), Bouchbinder, Look- 
man, and I proposed a statistical-thermodynamic frame¬ 
work for studying dislocation-mediated plasticity in poly¬ 
crystalline solids. Our purpose was to replace, or to find 
a firmer basis for, the largely phenomenological equa¬ 
tions of motion that have been used for many decades in 
attempts to describe the dynamic behavior of these ma¬ 
terials. We succeeded in computing the measured flow 
stress for Cu over fifteen decades of strain rate, and for 
temperatures between room temperature and about one 
third of the melting temperature, using just a few plausi¬ 
ble, physics-based parameters. We also computed stress- 
strain - i.e. strain-hardening - curves as functions of 
strain rate and temperature. The success of these efforts 
has given some credibility to our first-principles reformu¬ 
lation of dislocation theory. To test its limits of validity, 
we need next to look at specific phenomena and, more 
generally, to make closer contact with existing knowledge 
in this field. 

The LBL analysis is based on fundamental principles of 
statistical physics that, as pointed out in the introductory 
sections of that paper, are almost completely ignored in 
conventional dislocation theories. It seems obvious, how¬ 
ever, that complex behaviors involving chaotic motions of 
very large numbers of entangled dislocations can be de¬ 
scribed only in statistical terms. The equations of motion 
for these systems must necessarily be consistent with the 
laws of thermodynamics, especially the requirements of 
energy conservation and non-decreasing entropy. Thus, 
these equations must describe the flow of energy and en¬ 
tropy through driven, nonequilibrium systems. In LBL, 
they are expressed in terms of two dynamical variables: 
the density of dislocations p, and an effective tempera¬ 
ture X that describes the system’s state of configurational 
disorder. These two variables determine responses to ex¬ 
ternal driving forces, and carry the memory of earlier 
deformations. 

At the next level of specificity, the LBL analysis as¬ 
sumes that plastic deformation in polycrystalline solids 
is caused by the motions of line defects, i.e. disloca¬ 
tions, and that these motions are determined primarily 
by the rates at which dislocation lines become unpinned 
from each other under the influence of elastic stresses and 


thermal fluctuations. LBL argued that pinning times are 
generally very much longer than the times taken by dis¬ 
locations to move from one pinning site to another, and 
therefore neglected viscous forces affecting the motions 
of unpinned dislocations. Moreover, LBL paid no atten¬ 
tion to details such as the distinction between “edge” and 
“screw” dislocations, or “cross-slip,” or “stacking faults,” 
or the like. Instead, it was assumed that such details 
would become relevant if and when it was necessary to 
make first-principles estimates of quantities, such as rate 
factors, that appear in the equations of motion for p and 
X- The latter details do not appear in this paper either. 
However, “grain boundaries” appear in important ways. 

In what follows, I use the LBL statistical- 
thermodynamic framework to look at two long-standing 
puzzles in solid mechanics, neither of which has proved 
amenable to the conventional phenomenological anal¬ 
yses. First, I discuss the abrupt upturn in stress that 
has been observed in Cu for small values of the strain 
but at large values of the strain rate.0 Second, I 
discuss what I believe is a related observation that strain 
hardening curves for Cu sharpen abruptly when the 
grain size becomes sufficiently small. Q Understanding 
these phenomena illustrates how specific mechanisms 
can be incorporated into the LBL theory where nec¬ 
essary. It also illustrates how certain aspects of the 
thermodynamic theory - not fully developed in LBL - 
become relevant to practical applicatons. Therefore, in 
what follows, I start by reviewing the thermodynamic 
basis of the LBL analysis. 

II. THERMODYNAMIC EQUATIONS OF 
MOTION 

This thermodynamic analysis has been discussed in 
detail in earlier publications devoted primarily to the¬ 
ories of plasticity in amorphous materials. For exam¬ 
ple, see [3, Q. For both the amorphous and polycrys¬ 
talline cases, the analysis starts by dividing the system 
into configurational and kinetic-vibrational subsystems. 
The configurational degrees of freedom describe the rel¬ 
atively slow,i.e. infrequent, atomic rearrangements that 
are associated with irreversible plastic deformation; the 
kinetic-vibrational degrees of freedom describe the fast 
thermal and vibrational motions of the atoms. 
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For the polycrystalline case, it is useful to think of a 
slab of material lying in the plane of an applied shear 
stress. The dislocations oriented perpendicular to this 
plane are driven by the stress to move through a “forest” 
of dislocations lying primarily in the plane, thus produc¬ 
ing shear flow. Let the macroscopic area of this slab be 
A, and let its thickness be a characteristic dislocation 
length, say, L. 

The total internal energy of this system is 

Utotal = UciSc, p) + Ur{Sr). ( 2 . 1 ) 

Here, Uc{Sc,p) is the configurational energy of a poly¬ 
crystalline material containing dislocations, grain bound¬ 
aries, and other structural irregularities, p is the areal 
density of dislocations or, alternatively, the total length 
of dislocation lines per unit volume. SciUc,p) is the 
entropy of the configurational subsystem computed by 
counting the number of configurations at fixed values of 
Uc and p. Ur{Sr) is the kinetic-vibrational energy of 
this system, whose entropy is Sr. The kinetic-vibrational 
subsystem serves as a thermal reservoir. Its temperature, 
ksT = 9 = dlJR/dSR, is assumed to be fixed. The effec¬ 
tive temperature of the configurational subsystem, 

plays a central role in this analysis. 

Assume that we can write 

Uc{Sc,p) = Uo{p) + UiiSi), (2.3) 

and, correspondingly, 

Sc{Uc,p) = Sq{p) + Si{Ui), (2.4) 

where, Ui and are, respectively, the energy and the 
entropy of all the configurational degrees of freedom other 
than those associated with the dislocations. 

Uo{p) is the dislocation energy which, for present pur¬ 
poses, I write simply in the form 

Uo{p)=ApeD-, eD=L'-fD, (2.5) 

where eu is the energy per dislocation and is the dis¬ 
location energy per unit length. An implicit assumption 
here is that dislocations of opposite signs are present in 
equal numbers, thus minimizing the elastic energy. I also 
assume that the residual elastic energy, of the order of 
p ln(/9), is included approximately in eu, but is too slowly 
varying a function of p to be needed for present purposes. 
It may eventually be needed for studying spatially vary¬ 
ing dislocation patterns. Finally, in Ea. (l2.4|) . the entropy 
of the dislocations, S'o(p), can be written approximately 
in the form 

So{p) Ki — A p p) + A p for a^p<Cl, (2.6) 
where a is a length of the order of the atomic spacing. 


The usual thermodynamic analysis for this system goes 
as follows. The first law is: 

utotal = VaiP^ = Uc + UR 


= xSc+{^^ ^ P + eSR, (2.7) 

where V = LA is the volume, a is the shear stress, and 
is the plastic shear rate. (Variations of the reversible 
elastic energy cancel out of this equation. For example, 
see [^.) Use Eg. (12.71) to evaluate Sc^ and write the sec¬ 
ond law in the form 

= - W + f I - > 0, (2.8) 

A \ Xj 

where 

is the difference between the power delivered to the sys¬ 
tem and the rate at which energy is stored in the form of 
dislocations. 

Equation (12.81) is the sum of independent inequalities 
that, according to an argument originally due to Coleman 
and Noll, must be satified separately. Non-negativity of 
the term proportional to Sr implies that the heat flux 
Q, defined here to be positive when heat is flowing (as 
expected) from the conhgurational subsystem into the 
thermal reservoir, is 

Q = esR = K.{x-e), ( 2 . 10 ) 


where /C is a non-negative thermal transport coefficient. 

For present purposes, assume that the mechanical 
power, V a is always positive. Therefore, the remain¬ 
ing inequality is 

Use Eqs. (EH) and dH to write Uc = Uo + Ui{Sc — Sq), 
so that 

fdUc\ _ dSt) dF 

\ dp ) Sc dp ^ dp dp' 

F{p) = Uo{p)-xSo{p). (2.12) 

Equation (I2.II|) is satisfied by writing an equation of mo¬ 
tion for p in the form 




( 2 . 13 ) 


where A4 is a non-negative rate factor. Next note that 
Eas. (|2.5p and (12.6p imply that dF/dp = 0 when p = 
Pss (x), where 


Pssix) 


g-eo/x^ 


(2.14) 
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It is simplest to rewrite Eg. (12.131) in the linearized form 


p = M 



P 

Pssix). ' 


(2.15) 


The factor M must be proportional to the power per unit 
volume a , which is the only scalar rate in the problem. 
It has the dimensions of energy per unit volume per unit 
time. The left-hand side of Eg. (12.151) has the dimensions 
of length (of dislocations) per unit volume per unit time. 
Thus, writing this equation in the form 


p = Kp 


cr e 


■pl 


ID 



p 

Pssix). 


(2.16) 


is dimensionally correct and identifies the dimensionless 
factor Kp as the fraction of the input power that is con¬ 
verted into dislocations. The second term on the right- 
hand side of Eg. (12.161) can then be interpreted as the rate 
at which dislocations are annihilated as required by the 
second law. 

Having derived an equation of motion for p, return now 
to Eo. (l2.7p and rewrite this first-law equation in a form 
suitable for deriving an equation of motion for y: 


xSc = VaeP^ 



p-Q. 

Sc 


(2.17) 


Use the decompositions in Eos. (12.31) and (12.4p to write 
the left-hand side as 


. dS, dSo . dSo . 


which defines the effective specific heat Ce//. Next, make 
a similar expansion of the right-hand side of Eo. (l2.17p . 
and note that the term proportional to dSo/dp cancels 
out, leaving 

311 

VceffX = VaeP^-^p-Q. (2.19) 


forces, and not by the rate at which they occur. In LBL, 
we pointed out that the rapid rise in the stress seen at 
strain rates comparable to atomic vibration frequencies 
can be interpreted as a rapid rise in y when the system 
is being driven too fast for it to relax between rearrange¬ 
ment events. In the situations to be considered here, 
however, the strain rates are not so large, and therefore 
I consider only cases where Xss = Xo- 

To use this observation in Eg. (12.191) . note that y is 
comparable to the mesoscopically large energy cd, so 
that y ^ 0 and Q « /Cy in Eg. (12.101) . The require¬ 
ment that Xss = Xo tells us that fC = V a e^Vxo, so that, 
with Ea. (l2.5L the equation of motion for y becomes 


Ce// X = cr 



-IDP- 


( 2 . 20 ) 


In summary, the equations of motion for p and y, 
Eos. (12.161) and (I2.20p respectively, have been derived here 
using only basic principles of statistical thermodynamics 
and dimensional arguments. Equation (12.161) describes 
the flow of energy through the system of dislocations, as 
constrained by the second law of thermodynamics. Equa¬ 
tion (12.201) describes the flow of entropy; it is a restate¬ 
ment of the hrst law. 


III. DYNAMICS 

To complete this theory, we need relationships between 
the plastic strain rate the total strain rate e, and the 
shear stress cr. Start with the assumption that the elastic 
and plastic strain rates are simply additive quantities, so 
that 


(T =/r(e-e^'), (3.1) 

where p, is the shear modulus. The expression for the 
plastic strain rate is based on the Orowan relation: 

gP* = pbv, (3.2) 


In most STZ papers, my colleagues and I have neglected 
the analog of the term proportional to p on the right- 
hand side of Eg. (12.191) . because the energy content of 
STZ’s is negligible in comparison to that of all the other 
configurational degrees of freedom. Here, however, it ap¬ 
pears that a non-negligible fraction of the configurational 
energy may be stored in the dislocations. 

As argued, for example, in the steady-state value 
of the effective temperature, say Xss, must be a func¬ 
tion of only the strain rate e. (In steady state, the total 
strain rate e is the same as the plastic strain rate e^b) 
In other words, the steady state of disorder can depend 
only on the rate at which the system is being “stirred” 
by shearing. Moreover, we argued that Xss must go to a 
non-negative constant, say yo, in the limit of vanishing 
e. If e is slower than internal relaxation rates, then the 
steady state of disorder is determined only by the number 
of atomic rearrangements that are driven by the external 


where b is the magnitude of the Burgers vector and v is 
the average speed at which dislocations move across the 
system. That is. 


where £ = 1/y^ is the average spacing between dislo¬ 
cations, and l/Tp(tT) is the depinning rate. As in LBL, 
assume that depinning is a thermally activated process 
with a stress-dependent barrier of the form 

C/p(ct) = (3.4) 

Here, the the height of the unstressed barrier is defined 
to be ksTp, and the characteristic depinning stress ap 
is the Taylor stress. 


ctt = prh yfp, 


(3.5) 
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where ht is a reduced shear modulus of the order of 
/r/30. There is nothing sacrosanct about the exponential 
function in Eg. (13.411 . A linear approximation would be 
satisfactory, because the ratio a /ut turns out to vary by 
not much more than a factor of two in the experiments 
to be discussed here. 

With these ingredients, the depinning rate is 


1 

rp(cr) 


— /p(ct), 

To 


(3.6) 


where 


/p(cr) = exp 



(3.7) 


and To is an atomic time scale, set equal to 10“^^ seconds 
throughout this analysis. It is convenient to use tq to de¬ 
fine a dimensionless plastic strain rate, which, according 
to Ea. (l3.2p . is 


q{a, p) = To = y/p [/p(ct) - /p(-cr)]. (3.8) 

where p = b^p. Antisymmetry on the right-hand side of 
this equation is formally required; but the second term in 
the square brackets is negligible for positive stresses in or 
above the strain-hardening regime. Dropping that term, 
and setting ctt = can solve Ea. ()3.8ll explicitly 

for the stress as a function of the strain-rate: 


= In 


(JT 



In 



= iy{T,p,q). 


(3.9) 


It is useful to non-dimensionalize these equations, and 
to replace the time by the strain e as the independent 
variable, assuming a constant total strain rate e. Let 
a = 6 in the normalization of p given in Eg. (12.61) and in 
the definition of p in Eo. ()3.8ll . (As shown in LBL, this is 
equivalent to making a small change in the time scale tq.) 
Define x = x/^d- Also, let tq e = qo- Then, Eos. (l3.1L 
(12.161) . and (12.201) become, respectively, 


7 ^ ^ 

de [ go 


(3.10) 


^ ^ Kpb'^a q{a, p) / _ p \ ^ 
de xd qo V J ’ 


(3.11) 


dx _ O' P) (l ^ Y 

de CeffCD go [ Xo '’V e-i/xy 

(3.12) 

Finally, before looking at applications of this theory, 
note that we can simplify it by assuming that we are 
interested only in the plastic behavior during and after 
the onset of hardening. Because p is substantially larger 
than other stress scales in the problem, Eq. (I3.10p is a stiff 
differential equation that is accurately solved by writing 


q{a, p) = qo- Thus, Eg. (13.91) becomes an expression for 
the stress: 


<J = PT\/p ty{T,p,qo). 
Equations (13.111) and (13.121) become 


dp pt b'^ 
de 


Id 


'/p p(r,p,go) ( 1-r 


-i/x 


and 

de 


= n^y/p ty{T,p,qo) 


1 - ^ - Po ( 1 - 
Xo 


(3.13) 


(3.14) 


P 


=-i/x 


(3.15) 

where = pT/ceffen is a dimensionless number very 
roughly of the order of unity. 


IV. APPLICATIONS 


A. Onset of Strain Hardening 


Strain hardening is the transient approach to steady- 
state flow. The physical mechanisms that determine this 
transient are contained in the prefactors on the right- 
hand sides of Eos. (13.141) and (13.151) . especially in the di¬ 
mensionless conversion factor Kp. The steady-state be¬ 
havior is trivially independent of these prefactors. As 
shown in LBL, using p = - the steady-state solu¬ 

tion of Eg. (13.141) ~ in Ea. (l3.13p . produces temperature- 
dependent curves of stress versus strain rate that are in 
excellent agreement with experiment. That agreement 
is especially striking when xo becomes a function of the 
strain rate in the strong-shock regime. 

However, it is not entirely trivial to disentangle tran¬ 
sient from steady-state behaviors in evaluating Kp. To 
see this, look again at the onset of hardening as discussed 
in LBL. Near this onset, the dislocation density is rela¬ 
tively small, so that the second term in the parentheses 
on the right-hand sides of Eos. (13.Ill) and (13.141) is neg¬ 
ligible. It seems plausible that the stress in this regime 
is simply the bare Taylor stress, i.e. p multiplied by the 
strain required to move a dislocation a small fraction of 
an atomic spacing away from a pinning point, unmodi¬ 
fied by the thermal effects implicit in v{T, p,qo). If this 
conjecture is correct, then Eg. (13.Ill) becomes 

fAfr = Sflw (4^1) 

Id Id 

With no loss of generality, let xd = pb'^, where b' is 
a microscopic length, comparable to or perhaps smaller 
than the Burgers vector. Then 



1 ^ 1 /daT\ 

P V* Anset ~ P\ de Anset 


Kp p^ b^ 

2p^b'^ ’ 


(4.2) 


independent of p. According to Kocks and Mecking,0 
the hardening rate on the left-hand side of this equa¬ 
tion is often found experimentally to be about 0.05. If 
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firp ^ {b'/b) fi, i.e. if the length b' associated with the 
energy per unit length of a dislocation is the same as the 
displacement necessary to dislodge a dislocation from a 
pinning site in the formula for ctt, then Kp is about 0.1. 
In any case, if Hp is a constant of order unity or less, then 
the predicted value of the hardening rate in Ea. (j4.2l) is 
independent of strain rate and temperature, as observed. 

Equation (14.21) is not exactly what emerges when we 
use the full versions of Eos. (|3.131) and (|3.14|) to com¬ 
pute the onset rate. Instead, there appears an extra fac¬ 
tor u(T, p,qQ)^ multiplying Kp on the right-hand side of 
Eg. (14.21) . Thus, as in LBL, I conclude that Kp is propor¬ 
tional to k{T, p,qo)~^ in Eg. (13.141) . It seems to me that 
the physics leading to Eg. (14.21) is basically correct, as was 
the physics leading to the full equations in which Kp was 
an undetermined positive parameter. 


B. The Rate-Hardening Anomaly 

Materials scientists have been puzzled for decades by 
the sudden onset of rate hardening that is seen at strain 
rates of the order of 10"'’/sec, when stresses are mea¬ 
sured at small strains. The phenomenon is illustrated 
in FigU where the data points (taken from 0 and Q) 
show stresses measured in room-temperature Copper at 
four different strains, e = 0.05, 0.10, 0.15, and 0.20, as 
functions of strain rate. 

This behavior cannot be interpreted as an approach to 
some kind of singularity, as sometimes has been assumed 
in the literature. We know with fair certainty that the 
rate-hardening anomaly disappears if measured at larger 
strains, where the stresses must approach apparently un¬ 
remarkable steady-state values. Strain rates of 10"'^/sec 
are about four decades smaller than those that induce 
rapidly growing, disorder-generated hardening, as seen 
in the LBL analysis of the strong-shock regime. These 
two qualitatively different kinds of behavior should not 
be confused with each other. 

If the rate-hardening anomaly is a transient phe¬ 
nomenon, then its physical origin must be contained in 
the dimensionless prefactor Kp. It was assumed implic¬ 
itly, in the discussion following Eg. (j4.2l) . that Kp is a con¬ 
version factor that may incorporate a wide range of phys¬ 
ical mechanisms. One ingredient of Kp must be a rate at 
which dislocations are created, e.g., a density of Frank- 
Read sources or the like, multiplied by a nucleation rate 
per source, presumably proportional to the stress. To 
obtain a dimensionless conversion factor, multiply this 
rate by some time and volume associated with the driv¬ 
ing forces. The beauty of this approach is that we do not 
need yet to specify those details; we need only to know 
that Kp is a dimensionless number of the order of unity 
or less. 

The rate-hardening anomaly implies that Kp acquires 
a strain-rate dependence. With the contents of the next 
subsection in mind, I propose that this is a grain-size 
effect. Suppose that, whenever a grain of linear size d 



-4 -2 0 2 4 6 

Logio [Strain Rate (s”^)] 


FIG. 1: The rate-hardening anomaly as reported in and 
[^. The four curves, from bottom to top, show stresses 
as functions of strain rate for four different strains, e = 
0.05, 0.10, 0.15, and 0.20. The dashed curve at the top is the 
theoretical steady-state prediction. 


is sheared by an amount b/d, it emits a dislocation and 
then relaxes back toward its original shape. Thus, the 
emission frequency, multiplied by the areal density of 
grains, is proportional to e/db, in analogy to the source 
strength discussed in the preceding paragraph. This is 
an independent dislocation-creation mechanism that, to 
a first approximation, should be added linearly to the 
rate-independent term. 

This analysis, combined with that in Subsection IIV A1 
suggests that we write 


’^iT,p,qoy 



(4.3) 


where Rp is a constant of the order of unity, and qi is 
a dimensionless strain rate, possibly proportional to the 
grain size. By writing Eg. (14.31) in this way, qijrQ is ap¬ 
proximately the strain rate at which the upturn occurs. 
Then rewrite Eo. (l3.14p in the form 


— = K (i 

de ^^jz(r,p,go) V 



where 


(4.4) 


- b Pt 

Ki = Kp - 

ID 


(4.5) 


is a dimensionless number, again of the order of unity. 
Equation (13.151) is unchanged, but the factor Kp inside 
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the square brackets is now the function of j5 and qo given 
by Ea. (l4.3p . 

Figure [T] shows comparisons between experiment and 
theory for the rate-hardening anomaly. The dashed black 
line is the theoretical steady-state curve. The parame¬ 
ters used for plotting these graphs are almost the same 
as those used in LBL: Tp = 40800 if, T = 298 if, 
/iT = 1600 GPa, xo = 0.25, = 16, and ki = 3.1 

(the “universal” Kocks-Mecking value based on observed 
onset rates). This upturn analysis seems to be insensitive 
to the storage term in Eq. (13.151) ; therefore I have omitted 
it here by setting Kp = 0 in that equation. (That term 
will be nonzero in the next subsection.) The rate defined 
in Eg. (14.31) is qi/rQ = 4 x 10"^sec“^, with tq = 10“^^ sec¬ 
onds. The initial values of p and x used for integrating 
the differential equations are pi = 10“^ and Xi = 0.18, 
again the same as in LBL. 

Note that qi is the single new parameter that is needed 
to explain the rate-hardening anomaly for all four values 
of the strain shown here. The data for e = 0.15 is the 
least noisy of these data sets, and shows the cleanest 
upturn. It is this curve that appears most often in the 
literature. The agreement for the outlying cases, e = 0.05 
and 0.20 is somewhat less good, especially at the upturns 
where the data is most noisy. But the overall agreement 
between theory and experiment seems excellent. 


C. Effects of Grain Size 

As a second example of the statistical thermodynamics 
of dislocations, consider the effects of grain size in room- 
temperature Copper as observed by Meyers et al.Q That 
paper is largely devoted to measurements of grain-size ef¬ 
fects in phenomena such as dynamic recrystallization and 
strain localization. But it contains, in its Figs. 5 and 6, 
some data that is directly related to the present investi¬ 
gation. These figures show that the Kocks-Mecking rule 
of constant onset rate is violated at small grain sizes. The 
initial slopes of the stress-strain curves increase markedly 
with decreasing grain size. Moreover, at the smaller 
grain sizes, these curves exhibit narrow transitions be¬ 
tween rapid onsets at small strains and slower hardening 
at larger strains. The present theory can account for 
most of this behavior. 

Meyers et al. show stress-strain curves for materials 
with four different grain diameters: d = 9.5, 25, 117, 
and 315 pm. For each d, they show curves for strain 
rates 10“^ and 3 x 10^sec“^. To establish a base-line 
for the present analysis, start with the largest grain size, 
d = 315 pm, where there should be little, if any, anoma¬ 
lous strain-rate dependence. The theoretical results are 
shown in Fig|5]along with experimental points taken from 
0. The theoretical parameters are essentially the same 
as those used in Sec. IIVBI The only differences are that 
Ii/tq = 3 X lO^sec”^ (slightly smaller than before); and 
that = 12 and 8.5 for the smaller and larger strain 
rates respectively. Note that the onset rate is quantita¬ 



FIG. 2: Stress-strain curves for the largest grain size, d = 
315/rm . The strain rates are 10~® and 3 x lO^sec”^ for 
the bottom and top curves respectively. The data points are 
taken from [^j. 


tively as predicted by Kocks and Mecking, i.e. ki = 3.1, 
so that both curves have the same initial slope although 
their strain rates differ by more than six orders of mag¬ 
nitude. 

This situation changes dramatically when we look at 
the smaller grain sizes as shown in FiglH] Here, I 
have plotted stress-strain curves for three grain sizes, 
d = 9.5, 117, and 315 pm, all at the lower strain rate, 
do/To = 10“^, where the rate-dependent term qo/qi in 
Eg. (14.31) should be negligible. (I have omitted the inter¬ 
mediate curve for d = 25 pm for the sake of clarity.) The 
lower curve in this figure, for d = 315 pm, is the same as 
the lower curve in Fig|2] However, the upper two curves 
have substantially greater onset slopes, which means that 
the energy-conversion factor Hp or, equivalently, ki, must 
be a function of the grain size. The density of ordinary 
dislocation sources - not only the strain-rate driven ones 
discussed earlier - apparently increases with the density 
of grain boundaries. Accordingly, I have used the small 
strain-rate data in FigI3]to evaluate ki as a function of d. 
I then have used those values to compute the stress-strain 
curves at the higher strain rate, qo/'^o = 3 x 10^sec“^ 
shown in Fig 01 

The theoretical curves shown in Figs. 0] and 0] are 
uncertain in many respects, none of which seem to be 
fatal to the main concepts being tested here. In the first 
place, the data shown in Q is noisy, and my own ability 
to extract accurate points from it is limited. Secondly, 
there are too many things happening in this theory and, 
thus, too many parameters to evaluate. My strategy has 
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FIG. 3: Stress-strain curves at the strain rate 10“^ sec“^, for 
different grain sizes: d = 9.5, 117, and 315 fim, from top to 
bottom. The data points are taken from [^. 

been to start by fixing a few parameters based on prior 
experience with Cu data, e.g. in LBL and in the upturn 
analysis reported here, despite the fact that the materials 
used by Meyers et al are not exactly the same as those 
used elsewhere. Thus, I again use Tp = 40800 iT, T = 
298 K, fj,T = 1600 GPa, and xo = 0.25. The initial values 
of p and X are again pi = 10”^ and Xi — 0.18, except that 

1 have had to make small adjustments of Xi as noted 
below. 

The storage term, indicated by the factor Kp inside 
the square brackets in Ea. (l3.15L may be playing a role 
here. This is a thermodynamically predicted softening 
effect. The rate at which energy is stored in dislocations 
reduces the rate at which the entropy of disorder is in¬ 
creasing, thus reducing the rate at which the density of 
dislocations is increasing and, in turn, reducing the hard¬ 
ening rate. That mechanism seems to be playing a role 
here in reducing the slope of the stress-strain curves af¬ 
ter the initial rise at small grain sizes. I have tentatively 
accounted for it by letting kp = 0.2 in the formula for Kp 
in Eg. (14.31) for all grain sizes and strain rates. 

The values of the onset parameter determined at small 
strain rates, for d = 9.5, 25, 117, and 315/xm, are ki = 
9.5, 7.5, 5.0, and 3.1 respectively. The corresponding val¬ 
ues of Qi/tq, needed for the high strain-rate analysis, are 

2 X 10^, 6 X 10^, 7 X 10^, and 3 x 10^sec“^. Thus, as 
expected, both the ordinary conversion factor and the 
strain-rate induced one increase as the grain size de¬ 
creases. 

To fit the experimental data to the accuracy shown 
in the figures, I have had to make relatively small ad¬ 



FIG. 4: Stress-strain curves at the strain rate 3 x 10® sec“®, 
for different grain sizes: d = 9.5, 25, 117, and 315 pm, from 
top to bottom. The data points are taken from 0]. 


justments of other parameters. For the same increas¬ 
ing sequence of grain sizes listed in the preceding para¬ 
graph, I find = 18, 17, 17, 12 for the small strain 
rate, and 13, 11.5, 11.5, 8.5 for the large one. The cor¬ 
responding values of Xi are 0.165, 0.165, 0.17, 0.18, and 
0.16, 0.16, 0.16, 0.18. For d = 9.5 pm, the small strain- 
rate curve shown in Figl3]uses kp = 0.4 and xo = 0.253. 
So far as I can tell, these final adjustments serve only cos¬ 
metic purposes. Neither the experimental data nor the 
theoretical analysis are accurate enough for us to draw 
more definite conclusions from them. 


V. CLOSING REMARKS 

The analysis presented here contains several unconven¬ 
tional predictions that might usefully be checked exper¬ 
imentally. Specihcally, it predicts that the anomalous 
rate-hardening curves level off at higher strain rates, and 
that the anomaly disappears at larger strains. It also 
predicts that the upturn should occur at smaller strain 
rates for smaller grain sizes, and even speculates a simple 
power law dependence for this effect. In other words, it 
points toward extensions of earlier experiments. 

I close by repeating the plea from the conclusion of 
LBL. The theory presented there, and here, should be 
useful as a basis for studying dynamic, spatially hetero¬ 
geneous instabilities such as strain localization. If dislo¬ 
cation theory is to serve as a realistic tool for predicting 
the strength of materials, it must move in that direction. 






The thermodyamic STZ theory is developing successfully 
in that way as seen, for example, in its prediction of the 
fracture toughness of bulk metallic glasses. 0 I see no 
reason why the present dislocation theory cannot be used 
similarly. 
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